Entanglement of formation for a class of (2 Cg> d)-dimensional systems 
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Currently the entanglement of formation can be calculated analytically for mixed states in a 
(2 (g> 2)-dimensional Hilbert space. For states in higher dimensional Hilbert space a closed formula 
for quantifying entanglement does not exist. In this regard only entanglement bounds has been found 
for estimating it. In this work, we find an analytical expression for evaluating the entanglement of 
formation for bipartite (2 (g> d)-dimensional mixed states. 

PACS numbers: 03.67 Mn, 03.65 Yz 



Entanglement in multipartite quantum systems is an 
issue of fundamental interest in the fields of Quantum In- 
formation and Quantum Computation JTH5J. Although 
we know the entanglement of formation of a general 
mixed state of two qubits [Sj , for general bipartite mixed 
states in higher dimensional Hilbert space is unknown. A 
particular result is known for entanglement of formation 
of isotropic states in arbitrary dimensions [7] |5] . For gen- 
eral states the efforts have led to bounds for estimating 
the amount of entanglement [5HTT] . 

In this work we find an analytical expression for eval- 
uating the entanglement of formation for a class of 2 (g> d 
dimensional mixed states. Specifically we consider the 
calculation of EoF for 2 <Ei d bipartite mixed states ob- 
tained from tripartite 2®2®d pure states after tracing 
out one of the two-dimensional subsystem. Two physical 
models that meet such requirements are studied and the 
EoF is compared with a lower bound for EoF in higher 
dimensional systems [TU] , 

Our treatment appeals to recent developments in the 
study of quantum correlations for bipartite quantum sys- 
tems. The extension of classical mutual information to 
its quantum counterpart have shown that quantum cor- 
relations go beyond entanglement [T2]. The difference 
between quantum mutual information and classical cor- 
relation is defined as the quantum discord [13] . which 
measures the quantumness or Bayesian degree of a bi- 
partite state [Hj. A closed formula for quantum discord 
in a class of two-qubit states was achieved by Luo et 
al. [T5] and by Ali et al. [IB]. The study of such correla- 
tions has attracted much attention in last years [T7M22) . 
At this respect, Koashi and Winter found an identity for 
a tripartite (ABC) pure quantum state [23] , that relates 
the entanglement of formation (EoF) Eac of the reduced 
pair AC to the classical correlation (CC) J\~ B [23] [25J in 
the pair AB and to the von Neumann entropy of the re- 
duced state of A. Thus, the identity maps the problem 
of finding the EoF between A and C to the problem of 
finding the CC of the partition AB. This is schematically 
represented in Fig. [T] It is important to emphasize that 
the identity is not imposing restrictions on the dimen- 



sionality of the involved subsystems but on the purity of 
the tripartite state [25] . 
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FIG. 1. Scheme for the Koashi and Winter identity. Dou- 
ble line represents the classical correlation where an optimal 
measurement is carried out on system B. Single arrow line 
means entanglement of formation of the bipartite state pAC- 

Let us consider a tripartite pure state of systems A, 
B, and C, where A and B are qubits and C is a qudit 
(d-dimensional system). The correlations are distributed 
among the subsystems [23] as follows: 

Eac + Jab = Sa, (1) 

where Sa is the von Neumann entropy for the reduced 
state pa- The CC S24, 25] is defined by the optimization: 

JXb = max [S A - S(p AB \{IL x })] , (2) 

S (pab \{Kx}) — ~^2 x PxSa is the quantum conditional 
entropy of A averaged on all possible outcomes of the 
measurement {Tl x } performed on the subsystem B, p x — 
Tt{II x pab^x} is the probability with outcome x, and 
p x A = Trc{Tl x pAB^i x }/p x - Because of the nonsymmetric 
nature of the CC and by choosing three different ver- 
tices in Fig. [I] we find that there are six relations of the 
kind ([I]). Denoting the ordering ABC for Fig. [I] by even 
permutations we can obtain two additional equalities for 
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ordering CAB and BCA. On the other hand, by ex- 
changing B and C on Fig. [I] obtain an equality for the 
ordering ACB: Eab + JXc = ^A: which gives rise to 
equalities for BAC and CBA after permutations. Thus, 
according with equations and |2]) the EoF Eac is 
given by: 



Eac = min {S(pab\{Rx})} 



(3) 



Let us assume a general 2 ® 2 ® d tripartite pure state 
given by: 

I*) ABC = C 9g |g3) J 4B|V'l) + Cgel^ABl^) 

+ c eg |e.g)As|^3> + Ceelee)^!^}, (4) 

where \ipi) are arbitrary states of the <i-dimensional sys- 
tem C. By tracing out one qubit we are led with a 
(2 <g> (i)-dimensional mixed state. The relation ^ al- 
lows us to obtain the entanglement of formation for that 
bipartite mixed states, inasmuch as we can solve the 
optimization process to obtain the conditional entropy 
associated with the 2 ® 2 partition. The problem of 
optimizing 5(/j j 4b|{II i; } can be solved at least numeri- 
cally. However, there are some special classes of states 
Q where the 2 <g> 2 partition corresponds to a X state. 
This is the case for (-01 1-0 2 ) = (V'llV'a) — an d {i>4 l^fa) = 
(■04 1 ^3) = 0. For such classes an analytical recipe is 
available for solving the optimization process required to 
obtain Eac [23 HI- When these conditions hold, in the 
basis {\gg) = \G),\ge) = \l),\eg) = \2),\ee) = |3)} we 
lead with a state of the form: 



Tr c (\^)abc{^\abc) = 



( Poo p 03 \ 

pn P12 

P21 P22 

\ P30 p 33 ) 



(5) 



For this X state, the measurement on subsystem B leads 
to ensembles {po,po} and {pi,Pi} such that the quantum 
conditional entropy becomes: 



S(pab\{Rx}) = p S{po) +p 1 S(p 1 ), 



(6) 



where po, p\ are density matrices of the subsystem A 
after measuring the subsystem B, and pg, p\ are their 
respective probabilities. The corresponding eigenvalues 
are v±(p ) = (1 ± 0)/2 and io±(pi) = (1 ± 6>') /2 where 
6, 8' are defined in ref. [H]. The optimization has been 
reduced to a four parameter problem with four possible 
solutions which are extremes of the functional ^ |16j . 
We shall look for the set of parameters among these so- 
lutions which lead to the minimum value of the quantum 
conditional entropy. 

In what follows we consider two specific physical sce- 
narios where the state Q can be physically realized. Fol- 
lowing the recipes outlined in previous discussion we are 
able to calculate EoF analytically as a function of time 
for the mixed state of the (2 ® d) party. 



Our first example considers two atoms (A and B), 
each of them interacting resonantly with a single quan- 
tized mode of a cavity field (system C) in a Fock state. 
This physical situation is described through the two- 
atom Tavis-Cummings (TC) Hamiltonian: H — hg[(aA + 
cxb)clc + {°a + fT s) a c]: where Uj and crj are the ladder 
Pauli operators for the jth atom, a (a 1 ) is the annihila- 
tion (creation) operator for photons in cavity C and g 
is the coupling constant. Let us assume that the system 
is initially in the state |*(0)) = (a\gg) + (3\ee)) ab\ti) c ■ 
Since the TC Hamiltonian preserves the total number 
of excitation, the cavity mode will evolve within a 5- 
dimensional Hilbert space spanned by {|n — 2},|n — 
1), \n), \n + 1), \n + 2)} for n > 2. For n = 0, 1 the 
dimension will be 3 and 4 respectively. On the other 
hand, the atomic system will evolve within the subspace 
{\ gg ),\+),\ee)} with |+) = (\eg) + \ge))fy/2. By Solv- 
ing the Schrodinger equation, the system at time t is 
described by the state: 



|¥(t)> = a (t) \gg) \n + 2) + c 2 (t) |+) \n + 1) 
+ c 3 (t)\ee)\n) + c 4 (t)\gg)\n) 
+ c 5 {t)\+)\n-l)+c 6 (t)\ee)\n-2), 

where the probability amplitudes are 



(7) 



ci(t) = _ W(n+l)(n + 2) ^^-^^ (ga) 

c 2 (t) = - iP ^±} ^n(y/2(2n+3)gt), (8b) 
V2n + 3 

c 3(*)=^{ 1 - [l-cos(v^+3) S i)]}, (8c) 

c *(*) = «{ 1 - 2^1 [l-ooB(v/2(2n-l)st)] J, (8d) 



C5 * = - fi= a/2 2n-l )gt) , 
V2n— 1 



c6(t) = _a^g^I) J 1 _ cos( ^ 2(2n _ 1)5t) j (gf) 



The entanglement of formation £7ac between the atom A 
and the cavity mode C is given in terms of the quantum 
correlations embedded in the AB subsystem. The re- 
duced density matrix pab has the form of an X state ([5| . 
The matrix elements in that representation are p 00 = 
| Cl | 2 + |c 4 | 2 , p n = P22 = pa = P21 = (M 2 + |c 5 | 2 )/2, 
P33 = |c3 1 2 + I C6 1 2 , and p a3 = p^ = c 3 c%. Eac is the 
minimum of S(pab | {| n^}) in equation ^ and it 
is obtained analytically through the optimal eigenvalues 
{ v ±{0), w ±(0')} 01 the density matrices p and p\ [IB] . 

As an example, let us consider the case n — where 
the total system has dimension 2 ® 2 ® 3. In such case 
we found that there are two sets of parameters that 
can minimize Eq. ([6]), with corresponding eigenvalues: 
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FIG. 2. Entanglement of formation between the atom A and 
the cavity mode C for the initial state | \t/ (0) ) — (a\gg) + 
P\ee}) ab\0) c with a = /3 — l/v2. Blue line corresponds to 
analytical expression and black line to lower bound of EoF as 
a function of the dimensionless time r = \/6gt/2ir. 



FIG. 3. Entanglement of formation (blue line) and lower 
bound (black line) for a (2 ® 5)-dimensional state. The initial 
state is {a\gg) ab + /3\ee)AB) £S> |2)c with a = fi = 1/V2 and 
the dimensionless time r = \f\Agt/6-K. 



{v±( 



w±(9' 1 )} and {v±(0 2 ),w±(9 2 )}, where 



(a 2 



4) 2 



4(ac 3 + -|c 2 



\Ot + Ci 



^| C2 | 2 | 



62 = 

a 2 + ci + ±\c 2 \ 
At a time t the EoF is 



C 2 



ci + i| C2 i 2 



S AC = min{S'(^ ) ^) ! ^2 ) ^)}. 



(9) 
(10) 

(11) 



where we have written S(9,9') for the conditional en- 
tropy ^. For S , (6»i,6»i), po = p x = 1/2 and for 
S(9 2 ,9' 2 ), po = a 2 +4 + \c 2 \ 2 /2 and Pl = 1 - p - 
The dynamics of -E4C is shown in Fig. [2] where we re- 
stricted the dimensionless time r = V6gt/2n to one 
period. We observe that the entanglement is zero for 
t = 1/2 where the (2 <g) 3)-dimensional separable state 
given by p = \g)(g\ ® + /3 2 /9|e)<e| ® |0)(0|, with 

|0) = a|0) - 2/3V2/3|2). Although there is no entan- 
glement, this state is maintaining quantum correlations 
as can be checked by calculating the quantum discord 
which is different from zero except for a = (classical 
state) and a = 1 (factorized state). The amount of dis- 
cord in such case is called dissonance [14]. For r = 1 
the system is again in the factorized initial state. From 
Eqs. ^ we realize that the arising of entanglement is 
related with the population c 2 (t) of the symmetric state 
|+), e.L, entanglement will arise whenever 02(f) =/= 0. 

From Fig. [2] we see that the entanglement dynamics 
suffers abrupt changes. These happend for dimensionless 
times r « 0.40098 and r « 0.59902. These behavior 
arise because there are two solutions leading to equation 



(111. For r < 0.40098 the minimum conditional entropy 
is S(9 2 ,9' 2 ). For 0.40098 < r < 0.59902 the minimum 



conditional entropy is S(9i,9' 1 ). For r > 0.59902 the so- 
lution is again S(92,9' 2 ). These abrupt changes appear 
for all a different to zero. Evidence for the existence 
of abrupt changes in a lower bound for disentanglement 
dynamics in higher dimensional systems has been previ- 
ously reported 26 . 

By assuming the initial state |W(0)) = {a\gg) + 
fi\ee)) AB\n) c the atom- field system will evolve inside a 
(2 <g> 4)-dimensional Hilbert space when n = 1 and in- 
side a (2 ® 5)-dimensional Hilbert space when n > 2. 
Fig. [3] shows the evolution of the Eac for n = 2. In this 
case the symmetric state associated to the cavity state 
\n — 1) is also populated, i.e., entanglement arises when 
the probability amplitudes 02(f) and 05(f) are different 
from zero. We notice that at r = 1/2 the entanglement 
approaches to zero. At that dimensionless time the atom 
A and the field mode C are near the factorized state 
P = \9){9\ ® I^X^I wi th \4>) = p\2) + a\4). In other 
words, at t = 1/2 the state of the atomic subsystem 
has been approximately swapped to the field state. The 
minimum observed at r = 1 corresponds approximately 
to the initial state. On the other hand, for this initial 
state we found that there are three set of parameters that 
can minimize Eq. ([6]), that is, the EoF will corresponds 
to E AC = mm{S(9 1 ,9' 1 ), S(9 2 , 0' 2 ),S(6 3 ,6' 3 )} where the 
entropies can be obtained following the previous proce- 
dure. The dependence on three different entropies is the 
responsible for the abrupt changes showed in Fig. [3j 

Figs. |2] and [3] shows the lower bound [TU] used to com- 
pare with EoF. In both cases it follows monotonically the 
EoF, however it does not exhibit the abrupt changes. 

Different initial conditions for this system can be con- 
sidered whose evolution result in a X-type state for the 
atomic subsystem. In particular, initial entanglement 
between one atom and the cavity mode can be consid- 
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{v±(6 2 ),w ± (6' 2 )},with 



FIG. 4. Entanglement of formation between atom A and cav- 
ity mode (blue line) for the initial state |^(0)} = (a\gg) + 
/3|ee)) J 4s|0)c with a = 0.3. The red and green lines shows the 
evolution of the conditional entropies S(8i,8' 1 ) and S (02,02) 
respectively. Black line corresponds to the lower bound. 



ered. For instance 
|*(0)) =a\eg}\n + 



\y(0))=a\gg)\n + l)+p\ge)\n) and 
1)+P\ee)\n). 



A second physical scenario is a system consisting of 
two atoms A and B interacting with a common reservoir 
C. After tracing out the reservoir degrees of freedom the 
dynamics of the atomic subsystem is described by the 
master equation: p = (j/2)(2JpJ j; — jt J p — pj^ J), with 
J = a a + o\b • This model has been used to describe the 
dissipative dynamics of a collection of cold ions coupled 
through the center-of-mass mode [57]. However, by con- 
sidering explicitly the reservoir degrees of freedom and 
initially in the vacuum state the overall system will evolve 
in a (2 <g) 2 ® 3)-dimensional Hilbert space. Let us assume 
the initial state |^(0)) = (a\gg) + P\ee)) AB \Q)c, where 
|0) = Yik |0)fc is the reservoir the vacuum state. Inspect- 
ing the occupied Hilbert space throughout the evolution 
we found that the overall state is given by: 



|tf (t)> = a\gg} AB \0)c + c x (t)\ee) AB \U)c 

+ c 2 (t)\+) AB \l) c + c 3 (t)\gg) AB \2) c , (12) 



- a' - o-^t 



7 i — ci 



[/3 2 (2e 3 ^a + e 2 ^/3-/3 2 ) 5 



72 = 



47t_ 3e 27* /3 2 + ^2 )2] l/2 ; 

s 2 7 t / _ |2e 4 T*- 5e 2 T*/3 2 + 3^ 2 | 
3e*r*- 1' 02 = |2e 4 ^-3e 2 ^/3 2 +£ 2 | 



1 



with 
1/2 



In this case, the EoF is again given by Eq. 
Po = e- 4 T t (3e 2 ^*-l)/3 2 /2 for S(6l 2 ,6!!,) andp =pT 
for S^i, 6^). Fig. [4] shows the two relevant conditional 
entropies S(6i,6[) and S(82,0' 2 ), the EoF, and a lower 
bound. In this case the dynamics between atom A and 
the reservoir evolves in a (2 ® 3) Hilbert space. We ob- 
serve that there is only one abrupt change along the en- 
tanglement dynamics and that the atom-reservoir state 
converge to the factorized state \gg)(a\0) + VT — a 2 |2)) 
for time t >> 1/7. Notice that the reservoir converge 
to an coherent superposition when a / 0. The abrupt 
change holds at time for which the conditional entropies 
are equal and it depends on the initial probability am- 
plitudes. As in previous cases, the lower bound has a 
monotonic behavior with respect to the EoF and it does 
not exhibit the abrupt change. 

In summary, we have obtained an analytical recipe to 
evaluate the entanglement of formation for a class of 2®d 
bipartite mixed states. For this class of states the calcu- 
lation of EoF is mapped to an optimization process for 
calculating a conditional entropy in a 2® 2 partition. Pre- 
vious result for lower bounds shows to be in agreement 
with our findings. Finally, we have found abrupt changes 
in the entanglement dynamics as consequence of the op- 
timization process leading to the minimum value of the 
quantum conditional entropy. 

These results open the possibility for extending this 
treatment to evaluate EoF for arbitrary dimensions aided 
by such optimization for conditional entropies which can 
be realized numerically at least. 
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where \k) are collective states of the reservoir having fc 
excitations [28 . In this case the probability amplitudes 
become a(t) = /3e" 7t , 02(f) = fie^y/l - e'^r* and 
c 3 (f) = ^l-a2-cl(t)-4(t). 

This evolution takes to the atomic subsystem to be an 
X state where the matrix elements poo — c 2 , P03 — ctc\, 
P11 = P22 = P12 = P21 = cl/2 and p 33 = a 2 + c\. Here, 
we have also two relevant sets of parameters minimizing 
Eq. Q, leading to eigenvalues \v±(Q\) = w±(9' 1 )} and 
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